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Abstract 



We develop a quantum Smoluchowski equation in terms of a true probability dis- 
tribution function to describe quantum Brownian motion in configuration space in 
large friction limit at arbitrary temperature and derive the rate of barrier crossing 
and tunneling within an unified scheme. The present treatment is independent of 
| path integral formalism and is based on canonical quantization procedure. 
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^ ' 1 Introduction 

The theory of quantum Brownian motion is one of the major issues in physics 
and chemistry today. Its tremendous success in the treatment of various phe- 
nomena in quantum optics [1], quantum tunneling and coherence effects in 
condensed matter physics [2-4], activated processes in chemical physics [5-10] 
is now well-documented in the current literature. While the noise operator 
and density operator semigroup methods formed the core of development in 
quantum optics in sixties and seventies, real time path integrals attracted 
wide attention since eighties. However, these approaches are plagued with 
several difficulties. First, a search for quantum analogue of classical Fokker- 
Planck equation with a nonlinear potential, in general, leads us to equations 
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of higher (than two) derivatives of quasi-probability distribution functions 
[2,3,11]. These distribution functions often become negative or singular in the 
full quantum regime. Second, although large coupling constants and large cor- 
relation times are treated nonperturbatively formally in an exact manner by 
functional integrals their analytic evaluation often poses severe difficulties as 
emphasized recently by Stockburger and Grabert [12] and one has to consider 
special cases and approximations, e.g., WKB or semiclassical limit. Moreover 
because of well-known sign problem the numerical evaluation of the path inte- 
grals are very difficult [12]. Third, in some situations [13] the theory does not 
retain its validity as the temperature T — > 0. This implies that vacuum fluc- 
tuations due to heat bath are not correctly taken into account. The question 
is how to extend classical theory to quantum domain for large friction at arbi- 
trary temperature in terms of a true probabilistic description. Very recently we 
have addressed [10] this issue to develop a non-Markovian quantum Kramers' 
equation describing Brownian motion in a c-number phase space. The present 
analysis of overdamped limit is a brief offshoot of this development in config- 
uration space. Our aim here is (i) to develop a quantum analogue of classical 
Smoluchowski equation valid for arbitrary temperature and (ii) to obtain the 
rate of escape from a metastable state due to thermal activation and tunneling 
within an unified description. 



2 The quantum Langevin equation 



To start with we consider the standard system-heat bath model whose Hamil- 
tonian is given by 



h=- +V (x)+y: 



— H — K-j 
2 2 3 



(i) 



where X and P are the co-ordinate and momentum operators of the Brownian 
particle and the set {qj,pj} is the set of co-ordinate and momentum operators 
for the reservoir oscillators coupled linearly to the system through the coupling 
constant Kj and obeying the usual commutation relation [X,P] = ih and 
[QjiPj] = ifodij- The potential V(X) is due to the external force field for the 
Brownian particle. Eliminating the bath degrees of freedom in the usual way 
[1] we obtain the operator Langevin equation for the particle 

t 

X(t) + J dt' Pit - t') X(t') + V'{X) = F(t) , (2) 
o 
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where the noise operator F(t) and the memory kernel (3(t) are given by 

p(t) = ^[{^(0) - X(0)}kj cosujt +p j {0)K 1 / 2 smujt] and (3) 
j 

pit -f) = j2 K i cos (* - , ( 4 ) 

3 

respectively, with Kj = Very recently [10] we have shown that on the basis 
of a quantum mechanical average (...) over the bath modes with coherent 
states and the system with an arbitrary state Eq.(2) can be cast into the form 
of a generalized quantum Langevin equation in c-numbers, 

t 

x + J dt'(3(t - t')x(t') + V\x) = f\t) + Q(x, t) . (5) 
o 

In writing Eq.(5) we denote the quantum mechanical mean value of posi- 
tion operator (X) = x and Q(x,t) = V'((X)) — (V'(X)). c-number quantum 
Langevin force f(t) satisfies (f(t))s = and 

</(*)/(0>s = \ll (ooth ^) cosu.it - . (6) 



Here (. . )s corresponds to an ensemble average [10] over the quantum me- 
chanical mean values of the co-ordinates and momenta of the bath oscillators. 
Eq.(6) is the celebrated quantum fluctuation-dissipation relation [2]. Quan- 
tum noise fit) and quantum fluctuation term Qix,t) are due to heat bath 
and nonlinearity of the potential, respectively. For details we refer to [10]. 



3 Large friction limit and quantum Smoluchowski equation 



We now consider the diffusion of a particle in an external potential V(x) as 
described by Eq.(5). In the overdamped limit we drop the inertial term x 
and assume a Lorentzian density of modes of heat bath oscillators such that 
niuo)piuo) = (2/7r)[7/(l+cj 2 r c 2 )], which in the limit r c — > results in a damping 
kernel (3{t — t') = 7#(t — t') in Eq.(5). p(u) refers to the density function used 
in the continuum limit. Eq.(5) then assumes the form 

± + ±V^(x,t) = IQ (7) 
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where we have expressed the effective quantum potential V quant (x, t) as Vg Uant (x, t) 
= V'{x) — Q(x,t). Making use of the Stratonovich prescription the equivalent 
description of Eq.(7) in terms of true probability distribution p(x,t) is given 
by 

^ 18 [V^*,t)] + I>,.g. (8) 



dt 7 dx 



Here D qo is the quantum diffusion coefficient which can be obtained with the 
following definition [1] 

t+At t+At , 

2Dqo= At I dh I ^(f^f^s (9) 
t t ^ 



where the correlation function {f(t\)f{t2))s of the c-number quantum noise 
of the heat bath is given by Eq.(6). Making use of Eq.(6) in the continuum 
limit in Eq.(9) we obtain 

1 °° ( hui \ 

2Dqo = 2^fAt I ^"H^M^ y coth 2k^r) lu} • ^ 10 ^ 



Here I w is given by I w = // +A * dt\ // +A * dt 2 cosc<j(ti — t 2 ) which after explicit 
integration yields 



4 . ujAt 

up 



L = - 2 sin^. (11) 



The general expression for quantum diffusion coefficient in the overdamped 
limit is therefore given by 

2 1 °° 1 1 At 

D - = ^/l7^ M2sM + 1, ^ sin2 V^' (12) 

c 



where we have put the Lorentzian density of states for bath oscillators as stated 
earlier. n(u) is the average thermal photon number of the heat bath and is 
given by n(u) = l/[exp(7iu; /ksT) — 1]. The quantum diffusion coefficient can 
thus be calculated exactly by numerical integration over the bath frequencies 
u) for r c — > (the Markovian limit). The following two limiting situations are 
further relevant for the present analytic treatment: (i) In the high temperature 
limit ksT ^> fhj, the quantity hw[2n(u) + 1] in Eq.(12) reduces to k B T and 
the explicit integration results in D qo = k B T/^^ the usual Einstein's value in 
the static friction limit, (ii) In the vacuum limit on the other hand n(u) — > 
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and D qo may be obtained by considering the fact that the frequency depen- 
dence of the integrand in (12) [or (10)] except I u is flat (Markovian). This 
results in D qo ~ where Cj is an average bath frequency which may be 

approximately taken to be the linearized frequency of oscillation in the well co 
(i.e., uj ~ u ) since at T ~ the dynamics is dominated by population at the 
bottom of the well. For an arbitrary intermediate temperature we must rely, 
however, on the general expression (12) for the quantum diffusion coefficient. 

We now note that Eq. (8) contains quantum corrections to all orders. To make it 
more explicit we return to quantum mechanics of the system in the Heisenberg 
picture to express the operators X and P as 

X(t) = (X(t)) + 5X and P(t) = (P(t)} + 5P (13) 

(X(t)) and (P{t)) are the quantities signifying quantum averages and 5X 
and 5P are quantum corrections. By construction (5X) and (5P) are zero 
and they obey commutation relation [SX,SP] = ih. Using (13) in V'(X) 
and a Taylor expansion around (X) (= x) it is possible to express Q(x,t) as 
Q{x,t) = -E n >2^MV n+1 (x){SX n (t)) and V quant (x,t) as 

V quant (x,t) = V(x) + ]T ±-V n (x)(6X n (t)) (14) 

„>2 nl 

where V n (x) is the n-th derivative of the the classical potential which gets mod- 
ified by the quantum corrections. To solve quantum Smoluchowski equation 
(8) it is therefore necessary to calculate the corrections (5X 2 (t)), {8X 3 (t)) : 
etc. To this end we return to operator Eq.(2) and make use of the relation 
(13) and a Taylor expansion of the potential V(X) to derive the following 
coupled equations in the overdamped limit: 

1± + V > {x) + Y -V n+1 (x)(5X n (t)) = f(t) (15) 

n>2 U - 

(5Xn(t)) = --V"(x)(8X"(t)) (16) 



where n — 2, 3, .... A decisive advantage in the treatment of overdamped limit 
is noteworthy. The equations for quantum corrections {5X n (t)) are closed in 
contrast to those for the system without its surrounding [14,15]. If one is 
interested in the local dynamics around a point x (say, at the bottom or top 
of a potential well) the set of Eqs.(16) gets decoupled from (15) and one can 
obtain simple estimates of (SX n (t)) since V"(x ) can be treated as a constant 
in such cases. More generally however, Eq.(8) can be combined with Eq.(16) 
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to provide an extended phase space description in terms of a true probability 
function p(x, rj 2 , r] 3 ,...,t) 



dp(x,ri 2 ,r] 3 ,. . . ,t) 1 ^ d „ 1 d , 

dt = l^ 2 n Wn [V ^n P \ + - Yx [V{x) P \ 



7 dx 



-j V n+l(x)VnP 



n>2 



+ D 



qo 



d 2 p 
dx 2 



(17) 



where r\ n (= (5X n (t))) for n = 2, 3, . . . span the space of variables signifying 
quantum corrections around quantum mechanical mean position x. 

We now discuss the classical and vacuum limits of quantum Smoluchowski 
equation (8). As mentioned earlier in the classical limit, D qo reduces to Ein- 
stein's classical diffusion coefficient ksT/^y. At the same time Q(x,t) vanishes 
so that V guant (x, t) goes over to V(x) and one recovers the usual classical 
Smoluchowski equation. In the opposite limit as T — > 0, however, both quan- 
tum noise due to system and vacuum fluctuation originating from the heat 
bath make significant contribution. D qo in this limit assumes the form hwo/^'j. 
In this context we refer to a recent treatment [13] of large friction limit in quan- 
tum dissipative dynamics to point out that the latter theory does not retain 
its full validity as T — > since the quantum noise due to heat bath disappears 
in the vacuum limit. Another noteworthy feature of quantum Smoluchowski 
equation (8) [ or (17) ] is that unlike Wigner function based equations [11] it 
does not contain higher (than two) order derivatives of probability distribution 
function. The positive definiteness of this function is therefore ensured. 



4 Decay of a metastable state 



Based on the quantum Smoluchowski equation (8) we now consider the prob- 
lem of escape from a metastable state. Over the last two decades this has 
attracted a lot of attention in classical theory of thermally activated processes 
[5]. Since a quantum analysis naturally includes tunneling as an integral part 
of the problem we turn to this issue in search of an unified description of 
tunneling and thermal activation. 

To begin with we first take into consideration of the time-scale of free motion 
of the system, u (u 2 = V"). Since according to Eq.(16) quantum correction 
(SX n (t)) varies on a time-scale (nuj 2 / , ~f)~ 1 it is convenient to make an average 

of this fluctuation over a period and replace (5X n (t)) by an average (5X n ) = 
to Jq^ (5X n (t))dt. This implies quantum noise does not change significantly on 
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a time-scale during which the system relaxes so that V(x) does not depend on 
time explicitly. Eq. (8) can then be recast in the form of continuity equation 
dp(x, t)/dt = —dJ/ dx where J can be identified as a current which is given 
by 



. D dp 
J = -V'(x)p+- 7 f 
7 7 ox 



(18) 



Here we have put 



V(x) = V(x) + J2-V n (x)(5X") 



n>2 



71! 



(19) 



and D = 7-D go . We now consider a cubic potential of the form V(x) = 
— (A/3)x 3 +Bx 2 (A and B are positive constants), with a metastable minimum 
at x = which is separated from the true minimum by a finite potential bar- 
rier at x = x b . Linearizing the potential V(x) at the barrier top we calculate 
the stationary flux J (dp/dt = 0) around this point in the usual way, 



J 



co b y/D 



27T7 



p(0)exp(-E/D) exp 



-(l/2)V"\x b ){5X 2 ) b x b 
D 



(20) 



where u b refers to the frequency at the barrier top and E is the activation 
energy (E = V(x b )). Again considering an equilibrium distribution or zero 
current situation (J = 0) around x = we linearize the potential V(x) and 
derive from Eq.(18), the population in the well as given by 



n a 



UJQ 



p(0) 



(21) 



where u refers to the frequency in the well at x — 0. Based on the reason- 
ing given earlier, we can obtain an estimate of quantum correction in (20) as 

=x b = (^/2ct>b)[l + (2^/7)] which follows from the solution of Eq.(16) 



and a quantum average with minimum uncertainty state (SX 2 (0)) x=Xb = 
(h/2uj b ). Furthermore x b may be expressed as x b = \/2E/loq. With these sim- 
plifications the rate, k is given by the ratio J/n a : 



k = expl 

2tt7 FV 



-E/D) exp 



f2E/u G ){hA/2u b ) (l + ^ 



D 



(22) 



The expression (22) is the quantum rate of decay of a metastable state at any 
temperature. The quantum feature appears in two different ways. First, since 



7 



D may be treated as a quantum analogue of k B T in the Boltzmann exponential 
factor, which reduces to hwo/2 in the vacuum limit and to k B T in the thermal 
limit, the quantum nature of the heat bath is manifested through D. Second, 
the quantum nature of the system is pronounced through quantum correction 
terms which are entangled with the nonlinearity of the potential and vanishes 
in the classical limit. Thus in the classical limit (22) reduces to the familiar 
expression for the Kramers'-Smoluchowski rate of escape [16], 



2-7T7 



exp(-E/k B T) 



(23) 



In the vacuum limit, i.e., at T ~ 0, on the other hand we have 

'A(V2EM(h/2cu b ) (l + ^) 



k vac = exp(-2£; / ftwo) exp 

Z7T7 



(l/2)^o 



(24) 



The vacuum contribution corresponds to quantum tunneling at zero tempera- 
ture. It is important to recall that the problem of zero temperature tunneling 
was first successfully addressed by Caldeira and Leggett [2,3] in early eighties 
which subsequently initiated a major advancement in the field of macroscopic 
quantum tunneling [5]. The main results in the context of overdamped limit 
can be summarized as follows: First, the damping causes an exponentially 
strong suppression of tunneling rate. Second, for strong dissipation there is a 
large region where thermal and quantum fluctuations interplay so that the to- 
tal rate of decay increases because the thermal activation is supplemented by 
quantum tunneling at finite temperature. These features can be easily recov- 
ered in the present theory. The rate coefficient in the vacuum limit k vac shows 
that as the temperature approaches zero tunneling decay is almost exponen- 
tial. Furthermore, at finite temperature, the presence of factor (hA)/(Dub) 
in the exponential in the total rate coefficient in Eq.(22) implies an interplay 
of thermal fluctuation in D and quantum fluctuation due to (HA)/ub arising 
out of nonlinearity of the system. This leads to a effective reduction of the 
activation barrier from its classical value E resulting in an increase of the rate 
of decay. It is thus apparent that although the results are qualitatively com- 
parable the underlying physical picture in the two schemes (the path integral 
and the present one) are different. The tunneling decay at low temperature 
in the formal case is primarily determined by the nature of bounce solution 
evaluated by a semiclassical steepest decent method which depends on the 
effective path between the turning points. At the high temperature (i.e., in 
the classical regime) an unified multidimensional WKB approach is generally 
advocated to recover classical result [4,5]. The classical theory in the present 
formulation on the other hand is contained in k c i (Eq.(23)) which shows that 
the rate of thermal activation varies inversely with damping constant accord- 
ing to Kramers'-Smoluchowski limit. The quantum rate coefficient (22) thus 



8 



interpolates between high temperature thermal activation and zero tempera- 
ture tunneling. 



5 Conclusions 

In conclusion, we have considered the strong friction limit in quantum stochas- 
tic processes to develop a quantum Smoluchowski equation. The rate of decay 
of a metastable state which includes thermal activation and tunneling within 
a single scheme has been derived. Since classical Smoluchowski equation has 
a diverse range of applicability, e.g., in fluctuating barrier [17], in thermal 
ratchet and molecular motors [18-20], in tunnel diodes [21] and in the de- 
cay of periodically driven metastable potentials [22-24] we hope the present 
quantum analogue to keep its potential in similar issues. 

This work was supported by the Council of Scientific and Industrial Research 
(C.S.I.R.), Government of India, under grant No. 01/(1740)/02/EMR-II. 
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